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ABSTRACT 

Recent  developments  In  experimental  methods  for  mea¬ 
suring  the  characteristics  of  plane  compression  waves  are 
reviewed,  and  analytical  methods  for  Inferring  constitutive 
relations  from  measured  wave  profiles  are  discussed.  A  gene¬ 
ral  method,  requiring  only  pressure  or  particle  velocity 
measurements,  Is  proposed  that  Is  applicable  to  arbitrary 
waves  In  which  equilibrium  or  steady  state  may  not  obtain. 

A  summary  of  current  knowledge  of  constitutive  relations  ob¬ 
tained  from  plane  wave  experiments  Is  also  presented. 
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DETERMINATION  OF  CONSTITUTIVE  RELATIONS 
FROM  PLANE  WAVE  EXPERIMENTS 


I.  INTRODUCTION 

The  behavior  of  solids  under  Impact  loading  at  stresses  above  the 
linear  range  of  material  response  has  been  a  subject  of  study  since  about 
1872. ^  The  earliest  work  and  much  subsequent  work  has  been  done  with 
experimental  arrangements  that  produce  conditions  of  one-dime.islonal  stress, 
l.e.  rod  experiments. 

Since  the  Seco.id  World  War  Increasing  attention  has  been  given  to 
plane  waves,  In  which  the  strain  Is  accurately  one-dimensional.  Because  of 
this  strain  condition  deformation  of  the  material  takes  place  under  super¬ 
imposed  nydrostatlc  pressures  and,  at  the  higher  stress  levels,  at  elevated 
temperatures  resulting  largely  from  adiabatic  compression. 

The  most  distinctive  feature  of  plane  wave  compression,  however,  Is 
the  existence  of  shock  fronts.  In  these  fronts  the  strain  rates  can  be 
exceedingly  high  -  In  many  cases  exceeding  the  resolution  capabilities  of 
the  recording  Instrumentation.  Nonequilibrium  states  necessarily  exist  even 
In  a  steady-state  shock.  Thus,  a  relatively  complete  time-dependent  and 
temperature-dependent  constitutive  relation  Is  required  to  predict  plane-wave 
propagation;  by  the  same  token  a  wide  range  of  Information  about  material 
properties  can  be  obtained  from  plane-wave  experiments. 

The  complexities  of  shock  experiments  are  mitigated  considerably  by 
the  accuracy  with  which  the  one-dimensional  strain  boundary  condition  can  be 
produced,  In  rod  experiments  the  corresponding  one-dimensional  stress  condi¬ 
tion  Is  less  easy  to  verify  and  has  been  a  difficult  and  annoying  problem  In 
experiments  near  the  Impacted  end. 
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The  constitutive  relation  obtained  from  shock  experiments  Is,  however, 
that  pertinent  to  one-dimensional  strain,  and  it  is  probably  unreasonable  to 
expect  that  from  experiments  In  this  geometry  alone  one  can  Infer  a  completely 
general  constitutive  relation.  Nearly  all  experimental  techniques  measure 
only  the  stress  component  In  the  direction  of  propagation.  The  component  in 
the  direction  tangential  to  the  wavefront  can  usually  only  be  inferred. 

Where  hydrostatic  compression  data  are  available  comparison  with  shock  data 
permits  the  shear  stress  behind  the  shock  front  to  be  deduced.  Alternatively, 
measurement  of  the  states  obtaining  In  the  relief  of  stress  (via  a  rarefac¬ 
tion  wave)  from  a  shocked  state  provides  much  Information  about  the  yield 
stress  under  shock  conditions. 

To  a  considerable  degree  shock  waves  have  been  used  as  experimental 

tools  for  the  study  of  high  pressure  physical  phenomena  and  several  reviews 

have  appeared  In  recent  years  which  emphasize  that  aspect  of  shock  wave  phys- 
2-6 

Ics.  For  these  studies  stress  anisotropy  and  stress  relaxation  effects 
are  undesired  complications  and  are  frequently  Ignored.  Such  effects  are  not 
negligible  at  lower  pressures  or  where  wave  propagation  Is  the  principal 
Interest. 

In  this  article  we  focus  on  the  experimental  aspects  of  nonlinear 
plane  wave  propagation  In  solids.  Our  attention  Is  thus  directed  princi¬ 
pally  to  stresses  between  the  elastic  yield  point  and  stresses  one  or  two 
orders  of  magnitude  higher.  To  predict  wave  propagation  In  this  regime 
requires  a  time-dependent  constitutive  relation  for  one-dimensional  strain. 

The  converse  of  this  problem,  namely  how  to  deduce  constitutive  relations 
from  observations  of  wave  propagation,  Is  the  central  theme  of  the  paper. 

An  earlier  review  of  this  aspect  of  shock  wave  physics  has  been  given  by 
Karnes. ^ 

The  remainder  of  this  section  Is  devoted  to  br  ef  description  of 
some  of  the  most  Important  features  of  the  shock  transition.  Section  II 
reviews  methods  for  producing  plane  waves,  and  describes  recent  developments 
In  recording  techniques  for  observing  wave  behavior. 
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In  Section  III  analytical  methods  for  extracting  Information  about 
constitutive  relations  from  measured  wave  profiles  are  discussed.  In  that 
section  I  propose  a  rigorous  method  for  the  reduction  of  pressure- time  data 
that  does  not  require  the  usual  assumptions  of  steady  state  or  equilibrium, 
Isentroplc  flow.*  Section  IV  presents  a  review  of  current  knowledge  of 
constitutive  relations  derived  from  plane-wave  experiments. 


A.  Some  Properties  of  Shock  Waves 
1.  Jump  Conditions 

By  "shock  front"  one  means  a  compressive  (usually)  wave 
front  which  Is  steady.  That  Is,  a  shock  front  transforms  the  mechanical 
and  thermodynamic  state  of  the  material  to  one  of  higher  stress,  density, 
energy  and  mass  velocity.  The  Initial  and  final  states  are  equilibrium 
states  but  the  transition  region,  or  shock  front,  necessarily  Involves  non¬ 
equilibrium  states.  Because  the  transition  region  Is  steady,  no  mass, 
momentum,  or  energy  accumulates  within  It  and  conservation  laws  can  be 
applied  to  relate  the  Initial  and  final  states.  These  conservation  laws, 

O 

called  the  Ranklne-Hugonlot  jump  conditions,  can  be  written: 


VVo  ■ 

1  -  [(u,-u0)/(U-u0)] 

(1) 

P1  ■  Po 

■  P0<u-uo)(Vuo) 

(2) 

II 

o 

LxJ 

1 

C(P,+P0)/2](V0  -  V,) 

(3) 

In  these  equations  V  (=  p"^)  Is  specific  volume,  P  is  the 
component  of  normal  stress  In  the  direction  of  shock  propagation  (not  neces¬ 
sarily  an  equilibrium  stress),  E  Is  specific  Internal  energy,  u  is  mass 
velocity,  and  U  Is  the  shock  front  velocity.  Subscripts  o  refer  to  the  state 

*While  this  paper  was  in  preparation  I  learned  that  Mr.  Roger  Williams  had 

37 

made  the  same  discovery  Independently. 
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ahead  of  the  shock,  while  subscripts  1  refer  to  the  shocked  state. 

It  is  readily  shown  that  these  equations  imply  that  the 
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entropy  of  the  shocked  state  is  higher  than  that  of  the  initial  state.  ’ 

The  locus  of  equilibrium  states  given  by  Eq.  (3)  is,  for  a  given  material 
in  a  given  initial  state,  a  unique  function  called  the  Hugoniot,  or  R-H 
curve,  which  therefore  lies  above  the  isentrope  through  the  initial  state. 

A  complete  equation  of  state  for  a  material  can  be  characterized  by  a 
family  of  R-H  curves  ce:  >.ered  on  different  initial  states. 

Note  that  the  components  of  normal  stress  parallel  to  the 
shock  front  do  not  enter  the  equations  directly.  They  influence  the  shock 
only  through  the  equation  of  state,  which  properly  should  be  termed  the 
equilibrium  constitutive  relation  for  one-dimensional  compression. 

The  jump  cond.tions  apply  not  only  to  the  equilibrium  end 
states  but  throughout  the  transition  region  since  each  portion  of  the  front 
is  also  steady,  equations  (1)  and  (2)  can  be  combined  to  give: 

u-“o  *  vo  ■'('W/'VV  (4) 

Since  all  parts  of  the  wave  travel  with  the  same  velocity,  U  -  uQ,  with 
respect  to  the  undisturbed  material,  the  locus  of  P,  V,  E  states  in  the  tran¬ 
sition  must  lie  on  the  straight  line  joining  the  initial  and  end  states  in 
the  P-V  plane.  This  line  is  called  the  Rayleigh  line.  (Figure  1.) 

The  difference  between  the  Rayleigh  line  and  the  R-H  curve 
at  a  given  volume  is  approximately  the  non-equilibrium  stress  obtaining  in 
the  transition  region  and  is  primarily  responsible  for  the  entropy  produc¬ 
tion*  If  the  material  is  treated  as  a  viscous  fluid,  the  steady-state  shape 

of  the  shock  front  can  be  derived  by  relating  the  non-equilibrium  stress  to 
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the  stress  rate  or  strain  rate.  ’ 


2.  Stability  of  Shock  Waves 
★  ‘ 

The  thermal  part  of  this  difference  is  normally  negligible  in  solids. 
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In  solids  a  single  shock  front  Is  frequently  unstable  and 

a  compressive  wave  propagates  as  two  or  more  shock  fronts.  The  stability 

criterion  Is  derived  by  assuming  the  shock  to  consist  of  two  fronts  and  com' 

2 

paring  their  relative  speeds. 

Assume  that  the  first  shock,  travelling  with  velocity  Ug 
In  laboratory  coordinates,  transforms  the  state  to  Pg  -  P^,  Ug  -  u^.  Its 
velocity  with  respect  to  the  material  behind^  It  Is  then  Ug  -  Ug.  The  velo¬ 
city  of  the  second  shock  Is  In  laboratory  coordinates,  and  -  Ug  with 
respect  to  the  material  ahead  of  It.  The  final  state  Is  P^,  V^. 

Employing  the  jump  conditions,  Eqs.  (1)  and  (2),  these 
velocities  can  be  written 


UB  "  UB  *  VB  /<VPA)/(VA-VB) 
uc  •  UB  =  VB  ^PC"PB^VB"VC^ 

If  the  assumed  second  shock  travels  faster  than  the  first  a  single  front  Is 
stable.  Thus,  the  condition  for  stability  is: 

(pC_PB^VB"Vc)  >  (PB"PA^VA“V 


Graphically  this  means  that  the  Rayleigh  line  joining 
point  Pg,  Vg  with  Pq,  is  steeper  (more  negative)  than  that  joining  Pg,  Vg 
with  PA,  VA.  (Figure  2.) 

A  cusp  in  the  P-V  curve  as  indicated  at  point  B  in  Figure 
2  occurs  in  the  majority  of  solids  at  the  elastic  yield  point.  Because  the 
strain  is  one-dimensional,  shear  stresses  are  developed  by  a  plane  wave.  In 
an  elastic  solid  the  relation  is:^ 


★ 

Although  this  derivation  is  not  entirely  rigorous,  more  detailed  examination 
shows  the  result  to  be  correct. 
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t  ■  [(l-2v)/(2(l-v))]P 

where  t  Is  the  shear  stress  and  v  Is  Poisson's  ratio.  Clearly,  as  P  In¬ 
creases  so  does  t  until  th'  material  yields.  The  elastic  portion  of  the 
P-V  curve  has  slope, 


-(dP/dV)  -  K  ♦  (4/3 )u, 

where  K  Is  the  bulk  modulus  and  u  the  shear  modulus.  In  the  plastic  region 
above  the  yield  point  the  slope  Is  reduced  to  (neglecting  work  hardening): 

-(dP/dV)  -  K 


Accordingly  there  Is  a  range  of  shock  amplitudes  for 
which  the  stability  criterion  Is  not  satisfied.  Since  K  Increases  with 
pressure,  however,  the  unstable  region  Is  bounded  at  higher  pressures 
(Point  D  In  Fig.  2)  as  well  as  at  the  elastic  limit. 

Instability  can  also  result  from  phase  transformations. 

The  Isothermal  volume  discontinuity  of  a  first-order  phase  change  frequently 
corresponds  to  a  cusp  In  the  R-H  curve  and  the  consequent  separation  of  shock 
fronts  Is  very  suitable  for  experimentally  detecting  the  transition  and 
measuring  Its  pressure. 


3.  Reflections  at  Interfaces 

For  plane  waves  the  Interaction  with  a  boundary  of  dif¬ 
ferent  shock  impedance  Is  characterized  by  continuity  of  the  stress  normal 
to  the  boundary  and  the  mass  velocity.  For  this  reason  It  Is  convenient 
to  consider  the  relations  between  stress  and  particle  velocity  obtaining 
in  shock  transitions  and  in  rarefactions. 

The  shock  velocity  can  be  eliminated  from  Eqs.  (1)  and 


(2)  to  give 
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“l  ’  “o  -  ‘•'fVoNVV  <5> 

From  this  relation  a  family  of  curves  can  be  plotted  In  the  P-u  plane,  once 
an  equation  of  state  Is  given,  that  represents  the  locus  of  equilibrium  P, 
u  states  attainable  by  a  shock  transition  from  a  given  Initial  state.  Except 
for  the  end  states  the  transition  states  do  not  lie  on  this  curve  but  on  the 
straight  line  joining  the  end  states,  analogous  to  the  Rayleigh  line  In 
the  P-V  plane.  Those  curves  with  positive  slope  are  pertinent  to  forward¬ 
facing  shock  fronts  (l.e.  shock  fronts  travelling  In  the  ♦x  direction);  those 
with  negative  slope  pertain  to  backward -facing  waves. 

For  rarefaction  waves,  which  reduce  the  stress  and  accele¬ 
rate  the  material  In  the  direction  opposite  to  that  of  propagation,  the 

#  ft 

relation  between  stress  and  mass  velocity  Is  given  by  the  Rlemann  Integral  : 

U1  '  yo  “  4  /(<lP/oc)  (6) 

where  c  Is  the  local  sound  speed. 

This  relation  can  also  be  represented  In  the  P-u  plane  as 
a  family  of  curves,  and, as  for  shocks,  forward-facing  waves  are  described 
by  the  curves  with  positive  slope,  while  backward-facing  waves  are  described 
by  those  with  negative  slope. 

Mhere  the  effect  on  the  R-H  curve  of  the  entropy  change 
Inherent  In  the  shock  transition  Is  small  the  two  families  of  curves  are 
the  same  and  no  distinction  need  be  made.  All  transitions  from  a  given 
Initial  state  must  lie  on  one  of  the  two  curves  passing  through  that  state. 

For  example,  consider  the  reflection  of  a  forward-facing 
shock  In  material  A  at  an  Interface  with  material  B,  assumed  to  possess 
smaller  shock  Impedance  than  A.  (Figure  3.)  The  Initial  shock  Is  repre¬ 
sented  by  P1 ,u1  and  lies  on  the  P-u  curve  of  material  A  centered  on  (0,0). 

*Th1s  result  Is  a  consequence  of  assuming  constant  entropy  and  therefore 
does  not  hold  across  shock  fronts. 
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Reflection  of  the  shock  at  the  boundary  produces  a  backward-facing  rarefac¬ 
tion  in  A  and  a  forward-facing  shock  in  B.  The  common  state  at  the  interface 
must  li  on  the  intersection  of  the  appropriate  curves  centered  respectively 
on  (P^  j)  and  (0,0),  i.e.  the  final  state  is  (P^.u^). 

If  material  B  were  a  free  surface  the  reflected  rarefac¬ 
tion  would  have  carried  material  A  to  zero  pressure  and  the  free-surface 
velocity,  u^.  Note  that  if  the  p-u  curves  for  shocks  and  rarefactions  are 
the  same  the  free-surface  velocity  is  just  twice  the  particle  velocity  prior 
to  reflection.  This  result  is  frequently  used  to  infer  particle  velocities 
from  measured  free-surface  velocities  for  well-behaved  materials. 

The  P-u  plane  is  an  indispensable  tool  for  qualitative  or 
semi-quantitative  analysis  of  complex  wave  interactions  and  a  set  of  curves 
for  known  materials  is  to  be  found  in  virtually  all  experimental  laboratories. 


II.  EXPERIMENTAL  TECHNIQUES 

A.  Production  of  Plane  Stress  Waves 

The  principal  tool  for  producing  plane  stress  waves  for  the 
study  of  constitutive  relations  is  the  single  stage  compressed  gas  gun. 

1 2-14 

For  precisely  controlled  impacts  these  devices  are  at  present  unsurpassed. 

Existing  guns  vary  considerably  in  their  design;  nevertheless 
there  are  certain  common  features.  They  are  all  smooth  bore,  usually  having 
been  drilled  to  close  tolerances  from  a  solid  forging  or  casting.  They  use 
compressed  nitrogen  or  helium  as  the  driving  gas,  pressurized  up  to  about 
6000  psi.  Substantially  improved  performance  would  result  from  the  use  of 
hydrogen,  but  handling  and  safety  problems  have  discouraged  use  of  this  gas. 
Siegel  has  given  an  exhaustive  treatment  of  the  gas  dynamics  in  guns  of  this 


The  projectile  diameters  vary  from  2-1/2"  to  6";  the  barrel 
lengths  from  about  10'  to  100'.  The  velocities  achieved  vary  from  about 
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0.1  mm/ys  to  1.5  mm/ys.  Higher  velocities  are  exceedingly  difficult  to 
reach  In  a  single  stage  gun,  except  that  use  of  hydrogen  could  increase  the 
velocities  to  about  2  mm/ys.  Velocities  lower  than  0.1  mm/ys  can  be  achieved 
but  non-reproducible  frictional  losses  tend  to  make  the  projectile  velocities 
erratic  In  this  range.  Table  I  shows  some  pressures  produced  by  impact  at 
1.5  mm/ys. 


The  barrels  are  evacuated  ahead  of  the  projectile  to  prevent  a 
gas  cushion  from  distorting  the  wave  shape.  Hard  vacuums  are  evidently  not 
necessary;  Barker  was  able  to  detect  no  effect  at  residual  gas  pressures 
below  0.6  torr  in  the  Sandia  3  meter  gun.^  Impact  with  the  target  usually 
takes  place  an  inch  or  two  in  front  of  the  muzzle  to  provide  space  for  expan¬ 
sion  of  the  projectile  while  still  maintaining  a  maximum  degree  of  alignment 
of  projectile  and  target. 

The  tilt  angle  between  the  projectile  and  target  must  be  pre¬ 
cisely  controlled  if  true  plane  waves  are  to  be  produced,  particularly  at 
lower  impact  velocities.  The  angle  of  the  wavefront  with  respect  to  the 
target  surface  is  frequently  an  order  of  magnitude  or  more:  larger  than  the 
misorientation  of  the  projectile  with  the  target  because  of  the  large  dif¬ 
ference  between  wave  velocity  and  impact  velocity. 

Large  tilt  not  only  produces  two-dimensional  flow  but  reduces 
the  time-resolution  of  the  recording  instrumentation.  Thus  if  the  recording 
gauge  has  finite  dimensions  in  the  plane  parallel  to  the  impact  surface,  the 
time  resolution  achieved  may  be  controlled  by  the  time  required  for  the  wave 
front  to  sweep  across  the  gauge. 

Ir  practice  angular  misorientations  of  a  few  tenths  of  a  mi 1 1 i - 
radian  are  commonly  achieved.  A  tilt  of  this  magnitude  would  typically 
result  in  a  time  resolution,  for  a  gauge  whose  lateral  dimensions  are  1/4 
inch,  of  about  5  nanoseconds.  Since  this  is  comparable  to  the  resolution  of 
typical  fast  oscilloscopes  the  desirability  of  small  tilt  is  obvious. 
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The  recoil  of  guns  of  this  type  can  be  difficult  to  contend 

with  while  maintaining  the  precise  orientations  required.  Where  the  target 

is  mounted  rigidly  to  the  barrel,  vibrations  can  precede  the  projectile  and 

cause  anomalous  signals.  A  convenient  solution  incorporated  into  the  gun 

of  the  Shock  Dynamics  Laboratory  at  Washington  State  University  is  to  allow 

the  gun  to  move  freely  while  holding  the  target  rigidly  on  a  separate  sup- 
14 

port.  The  motion  of  the  gun  greatly  reduces  the  recoil  forces  and  does 
not  seem  to  adversely  affect  the  tilt  provided  the  barrel  is  maintained 
nearly  torque-free.  A  drawing  of  the  WSU  gas  gun  is  shown  in  Fig.  4.  This 
gun  has  a  4"  bore,  a  length  of  46',  and  operates  on  nitrogen  or  helium  to 
6000  psi.  It  is  typical  of  many  guns  of  this  type. 

Guns  that  use  gunpowder  as  the  propellant  are  in  limited  use. 

Although  they  can  be  shorter  for  a  given  velocity,  and  therefore  are  less 
expensive,  the  problems  of  cleanliness  and  high  recoil  forces  make  them 
somewhat  less  desirable  for  studies  of  constitutive  relations. 

17  18 

For  very  high  velocities,  two-stage,  light-gas  guns  can  be  used.  * 
The  velocities  of  these  guns  can  reach  8  mm/ys;  however,  the  projectile  diameters 
are  generally  less  than  two  inches.  This  is  a  severe  restriction  since  useful 
one-dimensional  information  is  obtained  from  the  target  only  before  the  effects 
of  free  lateral  surfaces  influence  the  wave  shape.  Accordingly  only  simple 
measurements  of  shock  velocity,  from  which  equation  of  state  points  are  deter¬ 
mined,  have  been  attempted  with  these  guns. 

High  explosives  have  been  used  to  produce  plane  waves  and  are 
useful  for  extending  the  pressure  range  available  from  single-stage  guns. 

These  systems  utilize  a  plane-wave  lens  to  produce  a  plane  detonation  wave 
at  the  target  interface.  In  some  cases  useful  information  has  been  obtained 
from  two-dimensional  steady  state  waves  produced  by  a  running  detonation.^ 

In  addition  to  producing  higher  pressures,  explosives  are  con¬ 
venient  in  that  precise  synchronization  with  recording  instruments  can  be 
achieved.  This  is  important  when  high  speed  cameras  are  used  because  they 
frequently  record  only  over  a  small  time  interval  and  cannot  be  triggered 
by  the  event  to  be  observed. 
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Disadvantages  to  explosives  are  that  the  plane  waves  produced 
are  less  precisely  controlled  than  is  the  case  with  guns.  Moreover,  a  con¬ 
tinuous  range  of  impact  pressures  is  very  difficult  to  achieve.  For  these 
reasons  explosives  are  most  frequently  used  for  equation  of  state  measurements 
and  for  the  study  of  other  high  pressure  physical  phenomena. 

Electron  beams  and  lasers  are  also  beginning  to  be  used  for  the 
19  20 

production  of  plane  waves.  ’  These  produce  stress  waves  through  thermo¬ 
elastic  coupling  when  radiant  energy  is  absorbed  in  the  target  material  in 
times  short  compared  to  the  transit  time  of  rarefaction  waves. 

Electron  beams  provide  much  higher  energy  densities;  with  elec- 

p 

tron  energies  in  the  range  0.2  to  5.0  mev  the  maximum  fluence  (energy  per  cm  ) 

2 

varies  between  about  35  and  300  cal /cm  .  The  corresponding  stress  pulses  pro¬ 
duced  depend  on  the  absorbing  material  but  are  typically  ten  to  several  hundred 

ns  in  duration,  and  attain  peak  stress  amplitudes  up  to  about  100  kbar.  The 

1 9 

area  over  which  the  flux  is  reasonably  uniform  is  a  few  square  centimeters. 

Lasers  operated  in  the  pulse  mode  give  comparable  pulse  dura¬ 
tions  and  irradiated  areas,  but  deliver  total  energies  one  to  two  orders  of 
magnitude  less  than  do  electron  beams. 

These  methods  are  most  useful  when  it  is  desired  to  observe 
stress  wave  propagation  under  high  temperature  conditions  or  when  very  short 
pulse  durations  are  desirable. 


B.  Measurement  Techniques 

Parameters  which  are  directly  pertinent  to  a  description  of  plane 
stress  wave  propagation  are  the  normal  stress  component  acting  across  a  plane 
perpendicular  to  the  direction  of  propagation,  the  density  or  strain,  the 
internal  energy,  the  mass  velocity,  and,  in  general,  two  phase  velocities. 

The  existence  of  two  phase  velocities  is  not  widely  recognized  and  indeed  for 
steady  shock  fronts  or  isentropic  rarefactions  the  two  velocities  are  the 
same.  In  general,  however,  they  are  distinct,  as  is  shown  in  Section  III. 


12 


Three  independent  relations  expressing  conservation  of  mass, 
momentum,  and  energy  can  be  written  relating  the  above  parameters;  Eqs.  (1)- 
(3)  are  those  relations  for  the  special  case  of  steady  flow.  Consequently, 
three  measurements  are  required  in  order  to  determine  the  state.  The  normal 
stress  components  parallel  to  the  wavefront  are  not  directly  involved.  They 
can  be  inferred  under  certain  conditions,  however,  by  comparing  the  compres¬ 
sion  and  relief  paths  in  the  P-V  plane  or  when  a  hydrostat  is  known 
independently. 

The  measurements  which  can  be  most  accurately  made  at  present 
include  the  normal  stress  jomponent  and  the  free-surface  velocity,  which  can 
be  related  in  many  cases  to  the  mass  velocity.  For  some  materials  mass  velo¬ 
city  can  be  measured  directly.  From  measurements  of  these  quantities  at 
two  or  more  locations  in  the  sample  one  or  both  of  the  phase  velocities  can 
be  obtained. 


A  variety  of  optical  techniques  using  high  speed  cameras  have 
been  invented  and  used  successfully  to  record  surface  velocities  (and  wave 
velocities).  These  include  argon  flash  gaps,  inclined  mirrors,  moving 
images,  and  optical  lever  techniques.  The  precision  of  these  methods  is 
usually  about  1%  in  wave  velocity  and  5%  in  free  surface  velocity.  The 
time  resolution  is  generally  about  10  ns. 

Electrical  methods  are  also  in  wide  use  and  include  electrical 
shorting  pins,  inclined  resistance  wires,  and  capacitor  microphones.  These 
also  measure  free-surface  velocities  with  a  precision  comparable  to  that  of 
optical  methods. 

A  transducer  technique  in  wide  use  at  present  is  the  quartz  gage, 

which  measures  the  pressure  at  the  interface  between  the  sample  and  a  disc  of 

21 

x-cut  quartz  plated  on  both  flat  surfaces.  At  stresses  up  to  about  25  kilo- 
bars  the  current  developed  by  the  piezoelectric  polarization  of  the  quartz  is 
proportional  to  the  stress  at  the  quartz-sample  interface.  Useful  results 
have  been  obtained  at  higher  stress  levels  but  50  kilobars  is  about  the  upper 
limit.  It  has  high  inherent  time  resolution  and  is  convenient  to  use. 
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All  of  the  above  techniques  suffer  from  a  common  drawback; 
namely,  the  observations  are  made  at  an  interface  at  which  there  is  usually 
an  impedance  mls-match.  Consequently,  in  order  to  infer  the  character  of 
the  undisturbed  wave  requires  an  impedance-matching  analysis  similar  in 
principle  to  that  mentioned  in  Section  I.  This  analysis  cannot  be  performed 
rigorously  without  knowledge  of  the  constitutive  relation  of  the  sample, 
although  useful  information  can  be  extracted  by  making  some  reasonable 
approximations  and  with  subsequent  iteration.  For  some  materials  such  as 
aluminum  the  quartz  gage  is  a  reasonably  good  impedance  match  and  the  analy¬ 
sis  is  less  sensitive  to  uncertainties  in  the  constitutive  relation  of  the 
sample. 

Detailed  descriptions  of  the  above  techniques  have  appeared  in 

several  review  articles  in  recent  years  and  the  reader  is  referred  to  those 
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for  a  more  thorough  presentation.  ’  ’  The  remainder  of  this  section  is 

devoted  to  more  recent  developments;  these  include  pi ezoresi stive  gauges, 
electromagnetic  velocity  gauges,  sapphire  gauges,  and  laser  interferometry. 

4 

1.  Piezoresistive  Gauges 

Manganin  wire  was  first  used  as  a  pressure  transducer  in 

24 

hydrostatic  apparatus  by  Bridgman  in  1911.  It  is  desirable  for  this  pur- 

I 

pose  because  it  exhibits  a  positive  pressure  coefficient  of  resistance  and 
at  the  same  time  a  very  small  temperature  coefficient. 

25  26 

In  1964,  Bernstein  and  Keough,  and  Fuller  and  Price, 
reported  experiments  in  which  a  fine  manganin  wire  was  imbedded  in  an  epoxy 
disc.  The  disc  was  used  much  as  is  a  quartz  gage;  it  was  placed  against  the 
free  surface  of  a  sample  and  the  change  in  resistance  monitored  as  the  pres- 
'■•re  pulse,  transmitted  into  the  epoxy  by  an  initial  pulse  in  the  sample, 
passed  over  the  wire.  (Figure  5.)  These  experiments  established  that  the 
fractional  change  in  resistance  is  linearly  proportional  to  pressure  up  to 
about  300  kbar. 

I 
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Numerous  dynamic  experiments  have  yielded  pressure  coef- 
27 

ficients  in  the  range 

(1/R)(AR/AP)  =  2.0  to  2.9  x  10"3/kbar 

_3 

The  statically  determined  value  is  2.6  x  10  /kbar. 

These  discrepancies  have  not  yet  been  fully  resolved;  the 
values  seem  to  depend  on  the  supplier  of  the  manganin  and/or  the  calibration 
technique.  For  this  reason  some  investigators  use  manganin  gauges  at  pre¬ 
sent  primarily  as  interpolation  gauges  between  pressures  established 

28  -3 

independently.  For  commercial  manganin  a  value  near  2.9  x  10  seems  to 

be  most  widely  observed;  little  or  no  temperature  dependence  has  been  ob¬ 
served  so  that  the  calibration  should  not  depend  on  the  material  in  which 
the  manganin  is  imbedded.  There  is  some  indication  that  there  is  a  hyster¬ 
esis  effect  so  that  the  coefficient  may  be  different  when  measuring  the 
compression  part  of  a  pulse  than  when  measuring  the  rarefaction  portion. 

It  is  uncertain  whether  this  effect  is  real,  however,  or  what  physical 
mechanisms  might  be  responsible. 

In  spite  of  these  difficulties  it  seems  reasonable  to 
expect  that,  as  development  proceeds,  a  reproducible  gauge  with  a  well- 
determined  coefficient  can  be  fabricated. 

Because  of  impedance  mismatches  between  the  sample  and  the 
insulating  material  in  which  the  gauge  is  imbedded,  the  gauge  used  in  this 
mode  has  the  same  limitations  mentioned  above  when  an  undisturbed  wave  pro¬ 
file  is  desired. 


More  recently,  experiments  have  been  performed  in  which 

27 

the  manganin  is  imbedded  directly  into  the  sample  material.  In  this  mode 
a  relatively  undisturbed  record  of  the  shape  of  the  pressure  pulse  is  ob¬ 
tained.  A  variety  of  thin  elements  have  been  developed  for  this  purpose. 
They  are  typically  0.001"  or  less  in  thickness  and  frequently  are  in  the 
shape  of  a  grid  in  order  to  increase  the  resistance  of  the  active  part  of 
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the  gauge  while  maintaining  small  lateral  dimensions.  It  is  desirable  of 
course  to  keep  the  thickness  as  small  as  possible  in  order  to  increase  the 
inherent  time  resolution,  which  is  dependent  on  the  reverberation  time 
through  the  thickness.  Small  lateral  dimensions  are  also  desirable  to  mini¬ 
mize  losses  in  time  resolution  due  to  tilt  of  the  wavefront  with  respect  to 
the  plane  of  the  gauge. 

An  example  of  one  such  grid  is  shown  in  Fig.  6.  This  gauge 
was  photoetched  from  manganin  sheet  relied  to  0.00075"  thickness.  The  lateral 
dimensions  of  the  active  element  are  1/8"  and  the  resistance  is  about  1.5 
ohms.  This  gauge  has  four  terminals  in  order  to  use  a  bridge  method  for 
recording  the  resistance  changes.  With  the  bridge  method  and  a  constant 
current  power  supply  the  effects  of  stretching  of  the  current  or  voltage  leads 
due  to  edge  effects  is  minimized. 

When  the  sample  to  be  investigated  is  an  insulator  a  gauge 

of  this  type  can  be  inserted  directly  into  the  sample  with  only  a  very  thin 

layer  of  cement  to  fill  the  voids  between  the  grid  elements.  If  the  sample 

is  a  conductor,  however,  thin  insulation  must  be  added  to  prevent  premature 

shorting.  Insulating  materials  such  as  mylar,  mica,  glass  and  Lucalox  have 

been  employed.  Of  these,  Lucalox  is  attractive  because  it  has  high  electrical 

breakdown  potential ,  and  it  has  high  shock  impedance  so  that  the  impedance 

27 

match  with  metals  is  improved  over,  say,  mylar.  Because  of  its  low  compres¬ 
sibility  the  change  in  capacity  between  the  element  and  the  sample  is  also 

minimized.  A  disadvantage  is  that  fabrication  of  thin  films  is  difficult; 

29 

plasma  sprayed  films  are  one  possible  solution. 

The  time  resolution  of  these  gauges  when  used  in  metal 
samples  is  somewhat  poorer  than  the  time  resolution  of  quartz  gauges,  for 
example,  principally  because  of  the  insulation  thickness.  If  the  total  gauge 
thickness  including  insulation  is  several  thousandths  of  an  inch  and  several 
reverberations  of  the  pressure  pulse  are  required  to  establish  equilibrium 
between  the  gauge  and  the  sample,  the  time  resolution  can  be  of  the  order  of 
30-50  nanoseconds. 


*  I 
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Recording  durations  of  gauges  of  this  type  are  normally 
several  microseconds,  and  for  laboratory  use  are  generally  larger  than  the 
times  for  which  one-dimensional  flow  can  be  maintained.  Careful  treatment 
of  the  leads  is  important  for  longer  recording  times  since  shearing  of 
these  is  the  usual  cause  of  premature  failure. 

Experiments  with  vapor  plated  manganin  grids  and  with 
other  materials  such  as  calcium,  lithium  and  ytterbium  show  considerable 

30 

promise  for  improving  the  low  pressure  sensitivity  of  gauges  of  this  type. 
Calcium,  for  example,  exhibits  a  pressure  coefficient  roughly  ten  times  as 
great  as  manganin  at  pressures  at  least  up  to  28  kbar.  However,  it  has  much 
higher  temperature  sensitivity. 

An  interesting  variation  on  this  method  has  been  reported 
31 

by  Bernstein,  et  al .  They  have  used  two  manganin  grids  oriented  respectively 
in  planes  parallel  and  perpendicular  to  the  direction  of  propagation.  Each 
grid  measures  the  stress  component  normal  to  its  plane  so  that  the  difference 
in  signals  is  a  measure  of  the  effective  yield  stress. 


2.  Electromagnetic  Velocity  Gauges 

When  the  sample  to  be  studied  is  an  insulator  an  electro¬ 
magnetic  technique  can  be  used  to  measure  mass  velocities  directly.  This 
technique  was  first  used  to  measure  detonation  parameters  in  explosives. 

Its  use  in  inert  solids  was  first  reported  by  Dremin  who  used  it  to  deter- 

32 

mine  the  behavior  of  glass  under  shock  loading.  Ainsworth  and  Sullivan 

33 

have  also  reported  extensive  measurements  on  rocks  up  to  30  kbar. 

The  idea  is  that  a  fine  wire  or  foil  imbedded  in  an  insu¬ 
lator  develops  a  potential  difference  when  it  moves  in  a  magnetic  field. 

Thus,  to  the  extent  that  the  foil  motion  is  the  same  as  that  of  the  insulator, 
a  measurement  of  the  voltage  across  the  wire  can  be  related  to  the  mass  velo¬ 
city  of  the  insulator. 
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The  experiment  is  shown  schematically  in  Fig.  7.  If  the 
magnetic  field  strength  is  B,  and  t  is  the  length  of  the  foil  perpendicular 
to  both  B  and  to  the  mass  velocity,  u,  Faraday's  law  of  induction  reduces  to: 

i  =  -(a/at)  /  B-dA  =  Bjiu 

where  £  is  the  emf.  If  we  insert  typical  values  into  this  expression  we 
find 


>  5 

c  =  500  gauss  x  1  cm  x  10  cm/sec  =  0.5  volt. 

A  field  strength  of  this  magnitude  is  readily  attainable  and  yields  an 
easily  measurable  signal  at  mass  velocities  of  interest. 

It  is  perhaps  surprising,  in  view  of  the  simplicity  of  the 
method,  that  it  has  not  been  more  widely  used.  It  cannot  be  used,  of  course, 
on  conducting  samples  because  eddy  currents  induced  into  the  sample  would 
distort  the  magnetic  field. 

One  might  suspect  that  the  polarization  field  induced 
into  an  insulator,  by  virtue  of  the  dielectric  constant  of  the  insulator, 
would  distort  the  electric  field  in  the  wire  or  the  magnetic  field  itself 
and  influence  the  measurement.  Detailed  examination  of  these  effects,  how¬ 
ever,  shows  them  to  be  negligible.^ 

Ainsworth  and  Sullivan  state  that  the  method  is  useful  up 
to  about  30  kbar.  However,  it  is  not  clear  why  such  a  limitation  need  be 
imposed;  Dremin  reports  measurements  up  to  about  400  kbar  in  glass. 

When  used  with  guns  or  explosives  precautions  must,  of 
course,  be  taken  to  prevent  the  motion  of  the  projectile  or  explosive  gases 
from  distorting  the  magnetic  field. 

The  precision  reported  by  Dremin  is  approximately  3%.  As 
with  manganin  gauges  thin  foils  are  desirable  to  achieve  high  time  resolution. 
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3.  Sapphire  Gauges 

Some  recent  work  has  been  devoted  to  developing  sapphire 
35 

as  a  transducer.  It  is  used  in  similar  fashion  to  the  quartz  gauge,  but 
depends  for  its  operation  on  the  change  in  capacitance  due  to  the  change 
in  dielectric  constant  and  to  the  reduced  electrode  separation  resulting 
from  shock  compression. 

The  current  developed  by  the  gauge  is  a  function  of  the 
mass  velocity  of  the  impacted  surface;  the  relation  is  linear  at  low  velo¬ 
cities  and  is  expressed  as: 

i  (t)  =  (VAU/fi.2)  [y  +  (e^U)]  u(t),  0  <  t  <  l/\i 

In  this  expression  i(t)  is  the  observed  current,  V  is  the  initial  applied 
voltage  (of  the  order  of  2  kilovolts),  A  is  area  of  the  disc,  i  is  the 
thickness  of  the  disc,  U  the  shock  velocity,  e..  ard  y  are  the  unstressed 
permittivity  and  the  rate  of  change  of  permittivity  with  mass  velocity. 

At  higher  impact  velocities  the  relation  becomes  nonlinear, 
but  can  be  readily  expressed  in  terms  of  measurable  constants  of  the  material. 

Sapphire  in  the  60°  orientation  seems  to  be  usable  at 
impact  stresses  up  to  100  kbar  in  the  sapphire.  Because  its  shock  impedance 
is  relatively  high  it  provides  a  reasonably  good  impedance-match  to  heavier 
metals,  such  as  iron.  The  principal  disadvantage  is  the  short  recording  time 
available  from  reasonable  crystal  thicknesses,  caused  by  the  high  shock  speed. 
This  time  is  typically  0.25  ysec. 

Other  materials,  such  as  ruby  and  Z-cut  quartz  have  also 

35 

been  examined  as  possible  gauges  of  this  type.  The  lower  yield  stress  of 
ruby,  however,  limits  its  usefulness  to  stresses  below  about  40  kbar.  Z-cut 
quartz  is  not  suitable  at  present  because  it  exhibits  internal  conduction 
and  noise. 
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4.  Laser  Interferometry 

Interferometric  measurements  provide  the  highest  time  and 
space  resolution  currently  attainable.  For  these  methods  a  laser  Is  not 
only  convenient  as  a  coherent  light  source  but  Is  necessary  In  order  to 
achieve  the  requisite  high  light  Intensities. 

Two  schemes  have  been  reported;  both  of  them  were  developed 
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by  Barker.  ’  The  first  of  these  uses  the  laser  In  a  conventional  Mlchel- 
son  Interferometer  arrangement  and  Is  shown  schematically  In  Fig.  8.  The 
portion  of  the  laser  beam  reflected  from  the  mirror  surface  of  the  specimen 
Is  ccmpared  with  that  reflected  from  a  stationary  mirror.  Interferometric 
fringes  are  thus  formed,  each  of  which  corresponds  to  a  displacement  of  the 
surface  of  one-half  wavelength  and  the  spatial  resolution  Is  therefore  of 
the  order  of  0.3  micron. 

The  laser  beam  Is  focused  on  the  specimen  mirror  surface 
In  order  to  minimize  the  effect  of  projectile  tilt.  This  surface  can  either 
be  a  polished  free  surface,  In  which  case  the  problems  of  relating  free- 
surface  velocity  to  mass  velocity  are  the  same  as  In  many  cf  the  techniques 
mentioned  above,  or  It  may  be  a  mirror  surface  plated  on  an  Internal  sur¬ 
face  of  a  transparent  specimen.  In  this  case  a  direct  measure  of  mass 
velocity  Is  obtained.  In  either  case  Impact  stresses  must  be  limited  to 
those  for  which  the  mirror  retains  Its  Integrity. 

When  measuring  the  displacement  of  an  Internal  surface  a 
correction  Is  required  for  the  change  of  Index  of  refraction  of  the  shocked 
"window".  The  relation  that  best  fits  current  data  Is  the  Gladstone-Dale 
formula:^® 

dp/p  *  dn/(n-l) 

where  p  Is  density  and  n  Is  the  Index  of  refraction. 

The  uncertainty  Introduced  by  lack  of  complete  Independent 
knowledge  of  the  density  In  the  shock  (which  Is  one  of  the  parameters  one 
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wisr.es  to  determine)  does  not  produce  serious  errors  because  the  density 
changes  involved  are  usually  small. 

The  principal  disadvantage  to  the  technique  described 
above  is  that  the  spatial  resolution  is  generally  too  high.  Consequently, 
for  mass  velocities  greater  than  about  0.2  mm/ys  the  fringe  frequency  exceeds 
the  capabilities  of  current  recording  systems  (approx.  600  MHz). 

By  means  of  a  clever  modification  of  the  above  technique 
the  space  and  time  resolution  can  be  adjusted  over  a  wide  range;  moreover, 
the  fringe  frequency  is  proportional  to  the  acceleration  of  the  mirror  rather 
than  to  its  velocity.  Each  fringe  then  corresponds  to  a  velocity  increment 
of  predetermined  magnitude. 

In  this  modification,  the  "velocity  interferometer  tech¬ 
nique,"  interference  fringes  are  formed  by  superposition  of  two  portions  of 
the  laser  beam  reflected  from  the  specimen  surface  at  different  times.  The 
earlier  signal  is  delayed  a  predetermined  amount  with  respect  to  the  later 
signal.  The  arrangement  is  shown  in  Fig.  9. 

The  operation  can  be  understood  by  referring  to  Fig.  10. 

If  the  time  through  the  delay  leg  is  t  =  t«  -  t,  and  the  distance  travelled 

by  the  mirror  surface  in  that  time  is  S,  then 

S  =  UT 

wnere  u  is  the  average  surface  velocity  over  the  interval  t.  The  signal 
reaching  the  photomultiplier  at  time  (tg  +  tc),  where  t£  is  constant,  is  thus 
composed  of  the  signal  reflected  at  time  t2  plus  that  reflected  at  time  t^ . 

If  the  velocity  of  the  surface  is  constant  in  time  the  separation  of  the  sur¬ 
faces,  S  =  x(t^)  -  x(t^),  is  constant  and  the  fringe  frequency  is  zero.  If 

the  surface  accelerates,  however,  fringes  will  appear  at  the  rate 


( x/2 ) (dn/dt )  =  (dS/dt) 


T(dtr/dt) 


21 


or,  since  n  is  zero  when  u  is  zero, 

u(t)  =  (a/2t)  n(t) 

where  n  is  the  number  of  fringes  counted. 

The  resolution  of  the  system  is  controlled  by  the  delay 
leg.  For  very  small  delays  the  number  of  fringes  per  velocity  change  is 
small  and  the  resolution  in  velocity  is  correspondingly  reduced.  With  a 
typical  delay  leg  of  10  nsec  and  a  wavelength  of  6328  K,  the  coefficient, 

du/dn  =  31 .64  m/sec/fringe. 

The  time  resolution,  on  the  other  hand,  is  equal  to  the  delay  time.  This 
can  be  seen  by  observing  that  the  technique  effectively  measures  the  sepa¬ 
ration  of  two  surfaces  displaced  in  time  by  x.  Consequently,  a  constant 
velocity,  for  example,  will  not  be  observed  as  constant  until  both  surfaces 
move  with  constant  velocity,  i.e.  until  the  specimen  surface  has  travelled 
with  constant  velocity  for  a  time  x. 

The  balance  to  be  struck  between  these  two  resolutions 
depends  on  the  experiment.  The  values  indicated  above,  however,  show  that 
reasonably  good  resolution  of  both  time  and  velocity  are  attainable. 

Although  laser  methods  are  somewhat  restricted  in  appli¬ 
cation,  the  high  time  resolutions  attainable  and  the  accuracy  with  which  they 
can  be  calibrated  makes  them  essential  tools  in  the  experimentalist's  reper¬ 
toire. 


III.  INTERPRETATION  OF  EXPERIMENTS 

As  indicated  above  current  experimental  techniques  in  general  provide 
measurements  of  pressure-time  or  velocity-time  histories  at  locations  fixed 
with  respect  to  the  material,  i.e.  in  Lagrangian  coordinates.  If  the  mea¬ 
surement  is  made  at  a  boundary  where  the  shock  impedance  changes,  as  in 


f 


22 

free-surface  measurements  or  quartz  gauge  measurements,  the  recorded  data 
are  characteristic  of  two  superimposed  waves  -  the  incident  wave  plus  the 
wave  reflected  from  the  boundary.  To  separate  the  effects  of  each  wave 
requires  knowledge  of  the  constitutive  relation  (in  general,  time-dependent) 
of  the  material.  Thus,  in  principle  the  derivation  of  a  constitutive  rela¬ 
tion  from  such  measurements  requires  prior  knowledge  of  the  constitutive 
relation  and  the  analysis  is  somewhat  circular. 

Valuable  information  can  nevertheless  be  obtained  through  a  series  of 
successive  approximations.  The  question  of  convergence  of  these  methods 
does  not  seem  to  have  been  treated  theoretically,  but  frequently  the  results 
are  not  sensitive  to  small  errors  in  the  assumptions.  Some  of  the  techniques 
do  not  suffer  from  this  uncertainty  -  notably  piezoresistive  pressure  trans¬ 
ducers  and  electromagnetic  velocity  gauges.  However,  these  techniques  have 
other  limitations  so  that  most  of  the  techniques  in  use  are  complementary. 

Once  stress-time  and/or  velocity-time  data  are  obtained  for  the  undis¬ 
turbed  wave  in  the  sample  the  question  arises  how  to  interpret  the  data  to 
derive  a  one-dimensional  strain  constitutive  relation.  If  the  compressive 
part  of  the  wave  is  steady  the  jump  conditions  (Eqs.  (l)-(3))  are  valid  and 
the  analysis  is  straightforward.  In  the  rarefaction  portion  of  the  wave 
one  usually  assumes  that  the  states,  while  not  steady,  are  nevertheless 
equilibrium  states  and  that  the  stress  is  a  function  of  the  density  only,  i.e. 
the  flow  is  assumed  isentropic.  The  jump  conditions  can  then  be  applied  incre¬ 
mentally  to  yield  a  stress-density  "isentropic"  relation  for  the  rarefaction 
part  of  the  wave.  The  Riemann  integral,  Eq.  (6),  is  just  the  integrated 
momentum  jump  condition  (Eq.  (2)). 

Unfortunately  these  assumptions  are  frequently  not  met  in  an  experi¬ 
ment  and  the  interpretation  is  accordingly  not  rigorous. 

An  alternative  means  of  interpretation  is  to  assume  a  constitutive 
relation  and  attempt  to  reproduce  the  experimental  observations  by  trial  and 
error  using  a  computer.  This  method  is  not  only  expensive  but  offers  no 
guarantee  of  uniqueness. 
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A.  Arbitrary  Waves 

It  is  not  difficult  to  show  what  kinds  of  measurements  are  neces¬ 
sary  and  how  they  should  be  interpreted  in  order  to  deduce  directly  the  stress 
density,  and  energy  states  obtaining  in  an  arbitrary  wave.  This  method 
requires  no  assumptions  about  the  form  of  the  constitutive  relation,  steadi¬ 
ness  of  the  wave,  or  extent  to  which  the  states  are  er  n'librium  states. 

The  equations  expressing  conservation  of  mass,  momentum  and 
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energy  are,  in  Lagrangian  coordinates: 


(aV/at)  -  VQ( au/ah )  =  0 

(7) 

(au/at)  +  VQ( aP/ah )  =  0 

(8) 

(aE/at)  +  (P/p0)(au/ah)  =  0 

(9) 

In  these  equations  h  is  the  initial  (undisturbed)  coordinate  of  a  particle, 
p  is  density  (pQ  is  initial  density),  P  is  the  stress  component  in  the  direc¬ 
tion  of  propagation,  u  is  mass  velocity,  and  E  is  specific  internal  energy. 


These  equations  are  generally  valid.  The  only  restriction  is 
that  heat  flow  has  been  assumed  negligible  in  Eq.  (9).  Thus,  within  that 
restriction,  the  equations  apply  equally  to  compressive  or  release  waves 
and  are  independent  of  any  assumption  about  the  constitutive  relation. 

We  define  two  phase  velocities  associated  with  the  wave  as  follows 
cu  =  (ah/at)u 
and 


Cp  •  (3h/3t)p 


since  h  and  t  are  the  independent  variables,  P  =  P(h,t)  and  u  =  u(h,t),  and, 
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dP  =  ( 3  P/3 1 )dt  +  (9P/9h)dh, 


also, 


du  =  ( 9u/9t)dt  +  (9u/9h)dh. 


Thus, 


_  _  9P/9t 
P  3  P/ 3  h 


00) 


.  _  9u/9t 

U  9u/9h 


01) 


Combining  Eqs.  (7)  and  (11): 


(pQ/p2)(9p/9t)  -  (l/cu)(9u/3t)  =  0 


or 


p 

(pq/p  )dp  -  (l/cu)du  =  0,  along  h  =  const. 
Combining  Eqs.  (8)  and  (10): 


(9U/9t)  -  (1/ PQCp) ( 3 P/ 9 1 )  =  0 


or 


dP 


»0Cp  du- 


along  h 


const. 


Similarly: 


dE  =  (P/pQcu)du,  along  h  =  const. 


Summarizing,  we  have 


dV  ■  -du/i>0cu 


(12) 


} 
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dP  =  p  c  du 
o  p 

dE  =  (P/p0cu)du 


(along  h  =  const)  (13) 

(14) 


These  equations  can  be  applied  incrementally  to  observed  wave  shapes  to 
relate  velocities  to  stress,  density,  and  energy.  Note  that  Eqs.  (12)  and 
(13)  are  equivalent  to  the  ordinary  jump  conditions  in  differential  form 
(in  Lagrangian  coordinates)  except  that  different  phase  velocities  must  be 
used.  Only  for  steady  state  waves,  discontinuous  shock  fronts  or  isentropic 
simple  waves  will  these  velocities  be  equal  in  general. 


The  mass  velocity  can  be  eliminated  from  Eqs.  (12)  and  (13)  to 
give  the  stress-density  relation: 

2  2 

dP/dp  =  (p  /p  )c  c  ,  along  h  =  const.  (15) 

Clearly,  this  reduces  to  the  acoustic  relation  whenever  the  phase  velocity 
c  =  c  =  c  .  Similarly,  we  can  write: 

r  u 

dE  =  (PdP/p  ^c  c  )  along  h  =  const.  (16) 

U  |J 


Implementation  of  this  method  in  experiments  requires,  at  a  mini¬ 
mum,  measurements  of  pressure-time  profiles  at  more  than  one  location  in  the 
specimen.  From  these  the  velocities  Cp  can  be  obtained  and  mass  velocity 
profiles  deduced  from  Eq.  (13).  Thence,  cu  values  can  be  obtained  and  the 
density  and  energy  calculated.  Simultaneous  particle-velocity  measurements 
would  be  valuable  and  might  enhance  the  precision  but  are  not  essential. 
Particle  velocity  measurements  alone,  without  pressure  measurements,  would 
require  an  iterative  procedure  to  deduce  P  and  Cp  and,  although  less  conveni¬ 
ent,  could  be  used. 

The  data  obtained  are  values  of  one-dimensional  stress,  density, 
and  internal  energy  for  each  part  of  the  wave  as  it  progresses  through  the 
material.  From  this  information  one  can  fit  various  models  to  obtain  a 
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general  relation  for  stress  as  a  function  of  strain,  energy,  strain-rate, 
history,  etc.  The  stress  components  parallel  to  the  wavefront  might  also 
be  inferred  from  a  variety  of  such  measurements  or  from  independent  hydro¬ 
static  data. 


The  relation  between  the  phase  velocities  can  be  derived  from 
cp  -  cu  =  ( ah/at)  -  (ah/at)u 

=  (ah/at)p  -  [(ah/at)p  +  (ah/ap)t(aP/at)u] 

(aP/at)u 
cp  (aP/at )h 


or 

(aP/at)u 

1  '  cu/cp  =  (aP/at)h  ; 


Similarly,  (17) 

( au/at) 

1  "  cP/cu  ”  (au/at)^ 

Consequently,  the  phase  velocities  will  be  equal  whenever  the  stress  is  a 
function  of  velocity  only,  or  when  the  wavefront  represents  a  discontinuity 
in  stress  or  mass  velocity. 

There  are  two  cases  when  stress  is  a  function  of  velocity  only  - 
steady  flow  and  isentropic  simple  waves.  This  is  shown  in  the  following  two 
sections.  We  can  conclude  immediately,  however,  that  whenever  P  =  P(u)  the 
lines  of  constant  phase  are  straight  lines  in  the  h,t  plane.  That  is,  cu  = 

cp  =  cp(p)  from  Eq*  (13)»  and  (3cp/3h)p  =  °* 


B.  Steady  State 


l 
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For  steady  state  shock  fronts  the  phase  velocity  c  is  the  same 

r 

for  all  parts  of  the  wave  and  Eq.  (13)  can  be  integrated  to  give  the  usual 
momentum  jump  condition 


P  - 


P0(u-u0)(u-uo) 


where  Cp  has  been  replaced  by  the  shock  velocity  with  respect  to  material 
coordinates,  U-uQ.  Further,  since  P  is  a  function  of  u  only,  we  have: 


_  du  faP/at 


dP I du/dP 


3P 

9h 


=  c. 


Thus,  Eq.  (12)  becomes 


dp  =  [p2/p0(U-uQ)]du 

Integrating  yields  the  continuity  jump  condition 

1  -  pQ/p  =  (u-u0)/(U-uQ) 

The  stress-volume  curve  is  seen  to  be  a  straight  line  joining  the  end  states 
(Rayleigh  line)  for,  integrating  Eq.  (15)  yields, 


U  -  uQ  =  VQ  /(P-P0)/(VQ-V)  =  const. 

Finally,  Eq.  (16)  yields  the  Rankine-Hugoniot  relation: 

E  -  E0  =  (1/2)(P+P0)(Vq-V) 

Note  that  this  derivation  does  not  require  that  the  states  be  equilibrium 
states  (except  insofar  as  steady  state  implies  equilibrium  end  states).  More¬ 
over,  if  the  shock  front  is  considered  to  be  a  discontinuity  then  c  =  c 

r  ^ 

and  the  jump  conditions  hold  even  if  the  flow  behind  the  wave  is  unsteady. 


C. 


Isentropic  Flow  in  Fluids 


r 
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The  development  of  non-linear  wave  propagation  theory  in  fluids 

O 

relies  heavily  on  the  method  of  characteristics.  This  method,  in  turn, 
depends  on  the  assumption  that  the  flow  is  everywhere  particle-isentropic 
and  that,  therefore,  all  states  are  equilibrium  states. 

If  to  Eqs.  (7)  and  (8)  we  add  the  relations: 

(3S/8t)h  =  0  (18) 


and 


P  =  P(p.S)  (19) 

where  S  is  entropy,  we  can  write: 

p 

(aP/8p)^  =  a  on  h  =  const., 

where  a  is  the  sound  speed  with  respect  to  spatial  (Eulerian)  coordinates. 
This  relation  allows  us  to  eliminate  the  derivative  in  p  from  Eq.  (7): 

(pQ/p2)(3p/9P)(9P/9t)  +  3U/3h  =  0 

or 

3P/3t  +  (p2a2/pQ)(3u/3h)  =  0  (20) 

Multiplying  Eq.  (8)  by  a, 

pa(3u/3t)  +  (pa/p0)(aP/3h)  =  0  (21) 

Adding  and  subtracting  these  gives 

[(3/3t)  ±  (pa/pQ) ( 3/3h)]P  ±  Pa[(3/3t)  ±  (pa/p0)(3/dh)]u  =  0 


The  original  equations  are  now  reduced  to  a  directional  derivative  and 
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separate  into  ordinary  differential  equations  along  characteristic  lines. 
Thus,  along  characteristic  curves  given  by 

dh/dt  =  ±  pa/pQ  (22) 

the  compatibility  relations  obtain: 

dp  +  padu  =  0  (23) 


These,  combined  with  Eq.  (17) 


(3S/3t)h  =  0 

can  be  solved  by  stepwise  integration  to  yield  a  solution  to  the  flow  for 
given  initial  and  boundary  conditions. 

Where  the  wave  propagates  into  a  constant  state  (simple  wave) 
entropy  is  the  same  throughout  the  flow  and  Eq.  (23)  can  be  integrated  to 
give  the  Riemann  integral: 


u  -  uQ  =  /  dP/pa.  (24) 

This  relation  is  frequently  used  in  the  analysis  of  experimental  data  for 
rarefaction  or  continuous  compression  waves.  By  measuring  the  propagation 
speed  of  pressure  or  mass  velocity  increments,  Eq.  (23)  can  be  integrated 
numerically  to  yield  pressure-velocity  and  pressure-density  relations.  Note 
that  this  method  is  a  special  case  of  the  general  method  outlined  in  Section 
III. A.  For  this  case,  as  for  steady  flow,  the  two  phase  velocities,  Cp  and 
cu»  are  equal  since  P  =  P(u). 

The  characteristics  method  breaks  down  at  shock  fronts  not  nec¬ 
essarily  because  the  shock  becomes  a  discontinuity  but  because  the  assumption 
of  an  equilibrium  equation  of  state  is  untenable. 
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D.  Relaxing  Elastic-Plastic  Solids 

For  elastic-plastic  relaxing  solids  the  assumption  of  particle- 
isentropic  flow  is  retained,  or,  more  accurately,  the  effect  of  entropy 
changes  on  the  pressure-density  relation  is  ignored.  Nonequilibrium  states 
are  allowed,  however,  by  assuming  that  the  plastic  strain-rates  may  be 
time-dependent,  and  that  the  stresses  depend  only  on  the  elastic  strains. 

The  plastic  strains  are  assumed  incompressible. 

The  relation  to  be  added  to  the  conservation  relations  is: 

(aP/n)  -  a2(aP/at)  =  -F(P,p)  (25) 

where  F  is  the  plastic  strain-rate  multiplied  by  twice  the  shear  modulus. 

The  sound  speed,  a,  is  here  the  elastic  longitudinal  velocity,  a  = 

/(a  +  2m)/p.  Combining  this  equation  with  Eqs.  (7)  and  (8)  in  the  same 
fashion  as  for  isentropic  flow  in  fluids  yields  the  characteristic  equa¬ 
tions 


c+: 

dP  +  Padu  = 

-Fdt, 

on  dh/dt  =  pa/pQ 

c. : 

dP  -  padu  = 

-Fdt, 

on  dh/dt  =  -pa/p 

V 

dP  -  a2dp  = 

-Fdt, 

on  dh/dt  =  0 

These  equations  have  been  applied  to  the  solution  of  plane  wave  propaaation 
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in  quartzite,  which  exhibits  pronounced  time-dependent  effects.  They  have 

also  been  applied  to  infer  the  plastic  strain-rate  at  the  elastic  yield 

40 

point  in  iron  from  measurements  of  the  decay  in  elastic  amplitude. 

At  the  elastic  shock  front,  when  an  elastic  precursor  wave  exists, 
the  jump  conditions  can  be  applied  to  relate  P  and  u  through  the  relation 
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Hence,  on  the  leading  c+  characteristic, 

dP/dt  =  -F/2 


Measurements  of  the  rate  of  change  of  peak  stress  thus  yield 
a  measure  of  the  plastic  strain  rate  at  the  elastic  front.  Since  no  time 
is  assumed  available  for  dislocation  multiplication  the  plastic  strain  rate 
in  turr  determines  the  velocity  of  propagation  of  mobile  dislocations. 

In  the  general  case,  use  of  this  theory  to  deduce  constitutive 
relations  from  experiments  requires  that  the  constitutive  relation  be  varied 
by  postulating  the  function,  F(P,p),  and  comparing  computed  wave  profiles 
with  observations  until  agreement  is  achieved.  A  possible  weakness  is  that 
"viscous"  stresses  are  not  permitted. 

The  relation  between  F  and  the  phase  velocities  c  and  c  is: 

p  u 

2 

dP  -  a  dp  =  -Fdt  (on  h  *  const) 


and 


dP  *  (P02/P2)cpcudp 


(on  h  *  const) 


Therefore, 


-F  =  [1  -  (aV/cpcuP02)]OP/3t)h 

or  if  the  speed  of  elastic  longitudinal  waves  is  expressed  in  Lagrangian 
coordinates: 


* 


a 


-  ap/pQ 


F  «  [(a*2/cpcu)  -  l](3P/at)h. 


and 
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This  relation  can  be  used,  within  the  context  of  the  assumptions 
implied  by  Eq.  (25),  to  determine  the  plastic  strain  rate  associated  with  a 
measured  pressure  profile  at  any  point  of  the  profile.  Note  that  vanishing 
of  the  plastic  strain-rate  implies  that  CpCu  =  a*c  and  that,  in  general, 
the  product  cpcu  is  less  than  a*^.  To  determine  F  from  experiments,  however, 
requires  independent  knowledge  of  a*.  It  is  not  clear  how  this  quantity  can 
be  unambiguously  measured  in  wave  experiments.  However,  acoustic  measure¬ 
ments  under  hydrostatic  pressure  could  provide  the  necessary  data. 


IV.  EXPERIMENTAL  RESULTS 

A.  Equations  of  State 

At  higher  shock  stresses,  substantially  exceeding  the  yield 
stress,  the  hydrostatic  component  of  the  stress  tensor  dominates  and  shock 
experiments  effectively  yield  a  hydrostatic  equation  of  state.  This  is 
the  Rankine-Hugoniot  relation  (Eq.  (3))  mentioned  earlier,  with  the  x-com- 
ponent  of  stress,  P,  replaced  by  the  mean  stress,  P.  At  these  pressures 
the  elastic  wave  is  either  small  enough  to  be  neglected  compared  to  the 
"plastic"  or  deformational  shock,  or  may  be  overtaken  by  the  plastic  shock, 
i.e.  a  single  shock  may  be  stable.  The  shock  rise-times  are  also  usually 
very  small  so  that  the  transition  is  effectively  steady  state  whether  or  not 
the  final  state  is  a  true  equilibrium  state.  Hence,  the  R-H  jump  conditions 
can  be  applied  to  relate  the  shocked  state  to  the  initial  state  and  the 
more  elaborate  analysis  of  Section  III. A.  i*  unnecessary. 

A  very  large  amount  of  data  has  been  generated  on  equations  of 
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state,  particularly  of  metals  and  rocks.  *  Where  comparisons  can  be 

44 

made  with  hydrostatic  data  the  agreement  is  satisfactory,  and,  in  fact, 

shock  data  are  used  to  calibrate  static  pressure  apparatus  in  the  higher 
45 

pressure  ranges. 

A  principal  uncertainty  at  the  present  time  is  the  temperature 

of  the  shocked  state,  which  can  only  be  calculated  with  independent  knowledge 

3  41 

of,  or  assumptions  about,  the  Gruneisen  ratio  and  the  specific  heat.  ’ 


i 
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Attempts  to  measure  temperature  of  the  shocked  state  by  means 

of  the  thermoelectric  effect  have  thus  far  been  unsuccessful.  Anomalously 

large  signals  have  been  observed  whose  physical  origin  is  not  understood. 

Some  success  has  been  achieved,  however,  in  measuring  free-surface  tempera- 
47  48 

tures.  ’  These  provide  at  least  a  consistency  check  on  some  of  the 
assumptions  employed  in  interpreting  equation  of  state  experiments. 

By  using  porous  samples  R-H  curves  centered  on  different  initial 

states  (V0»Eq)  can  be  determined  and  the  Gruneisen  coefficient  measured. 

Thouvenin  has  questioned  the  assumption  that  equilibrium  states  can  be  achieved 
50 

by  this  method,  but  later  work  by  Hofmann  et  al  indicates  that  they  are 
possible.^ 


B.  Constitutive  Relations 


At  shock  stresses  comparable  to  the  yield  stress,  the  effects  of 
stress  anisotropy  and  strain  rate  cannot  be  neglected.  Not  only  does  the 
elastic  wave  carry  a  significant  fraction  of  the  total  stress,  but  the  struc¬ 
ture  of  the  plastic  shock  and  the  rarefaction  wave  that  relieves  the  shocked 
state  are  more  complex.  Figure  11  shows  an  example  of  the  compressive  por¬ 
tion  of  the  wave  shape  in  single  crystal  Lithium  Fluoride.  It  would  clearly 
be  a  coarse  assumption  to  consider  these  waves  as  simply  two  discontinuities 
in  stress  separated  by  a  constant  region.  Much  more  detailed  analysis  of 
the  wave  structure  than  simple  application  of  the  jump  conditions  is  required. 

Where  strain-rate  effects  can  be  ignored, elastic-plastic  theory 

can  be  applied  to  predict  the  differences  between  the  one-dimensional  consti- 
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tutive  relation  and  the  hydrostatic  equation  of  state.  *  This  theory 
predicts  the  curves  shown  as  Fig.  12.  Above  the  yield  point  the  difference 
between  the  normal  stress  in  the  direction  of  propagation  and  the  hydrostatic 
pressure  is  just  (2/3) Y ,  where  Y  is  the  yield  stress  in  simple  compression. 
Work-hardening  can  be  Incorporated.  The  release  curve  representing  the  states 
through  which  a  rarefaction  wave  carries  the  material  is  also  offset  from  the 
hydrostat  by  (2/3) Y  once  the  yield  stress  in  the  reverse  direction  is  attained. 


r 
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11  53 

Several  attempts  have  been  made  to  verify  this  model.  *  Mea¬ 

surements  of  the  compressive  state  in  aluminum  indicate  reasonably  good 
agreement  as  indicated  in  Fig.  13.  The  rarefaction  portion  of  the  curve  for 
the  lower  range  of  shock  stresses  shows  poorer  agreement,  evidently  because 
of  Bauschinger  effect. 

The  Hugoniot  elastic  limit,  indicated  by  PQ  in  Fig.  12,  is  fre¬ 
quently  observed  to  be  time-dependent  and,  as  would  be  expected,  is  also 
dependent  on  the  history  of  the  specimen.  Nevertheless,  the  values  observed 

are  approximately  characteristic  of  the  material  and  a  useful  tabulation  of 
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measurements  to  date  has  been  given  by  Graham  and  Jones. 

At  shock  stresses  of  100  kbar  or  more  the  difference  between 
Hugoniot  and  hydrostat  is  difficult  to  resolve;  moreover,  good  hydrostatic 
data  do  not  exist.  Measurements  of  the  rate  of  decay  of  a  shock  pulse  of 
finite  width,  however,  permit  inferences  to  be  made  about  the  yield  stress 
and  shear  modulus  under  shock  conditions. 

The  first  experiments  to  observe  shock  decay  were  performed  by 
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Altshuler  and  Curran  on  aluminum.  Curran's  results  indicated  that  shock 
decay  is  more  rapid  than  would  be  predicted  by  the  clastic-plastic  model 
with  constant  yield  stress.  He  postulated  that  yield  stress  increases  with 
compression,  attaining  a  value  of  12  kbar  at  a  relative  volume,  V/VQ,  of 
0.86.  The  initial  yield  stress  was  0.5  kbar. 
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These  results  were  verified,  at  least  qualitatively,  by  Erkman. 

He  found  similar  behavior  in  copper  and  epoxy. 

Recently,  van  Thiel  and  Kusubov  have  measured  the  shape  of  the 
rarefaction  portion  of  an  initially  square  pulse  induced  in  aluminum  by  im¬ 
pact  with  a  thin  plate.28  They  used  manganin  gauges  to  interpolate  between 
the  peak  pressure  of  the  pulse  and  zero  pressure.  The  impact  pressure  (130 
kb)  was  higher  than  in  Barker's  similar  experiments.  Their  conclusion  is 
that  the  yield  stress  is  a  peculiar  function  of  pressure  in  the  rarefaction 
wave  and  approaches  values  as  high  as  28  kbar,  or  nearly  five  times  the  sta¬ 
tic  strength.  (Fig.  14) 
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Their  analysis  is  based  on  the  assumption  that  the  initial 
shocked  state  lies  on  the  hydrostat  rather  than  above  it  by  (2/3) Y  as 
predicted  by  elastic-plastic  theory.  Other  experiments  by  these  investi¬ 
gators,  involving  free-surface  velocity  measurements  at  an  impact  pressure 

of  313  kbar,  indicate  a  yield  stress  of  14  kbar  when  analyzed  according  to 
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conventional  elastic-plastic  theory.  The  two  sets  of  experimental  results 
are  compatible;  the  factor  of  two  discrepancy  is  the  result  of  different 
assumptions  about  the  shocked  state. 

Similar  experiments  have  been  conducted  in  magnesium  by  Fuller 
and  Price. ^  Their  results  are  also  compatible  with  the  assumption  that 
yield  stress  increases  under  pressure,  attaining  values  perhaps  two  to  five 
times  the  zero  pressure  value  (1  kb)  at  a  shock  stress  of  80  kb. 

A  limited  amount  of  work  comparing  shock  and  hydrostatic  com¬ 
pression  has  also  been  performed  on  covalent  and  ionic  solids.  The  first 
experiments  were  those  on  quartz,  which  appeared  to  show  that  the  stress 

difference,  (2/3) Y,  that  might  be  expected  between  the  R-H  curve  and  the 

fil  fi? 

hydrostat  above  the  yield  point  does  not  obtain.  ’  The  conclusion  was 

that  the  material  loses  rigidity  at  the  yield  point  and  that,  therefore, 

the  yield  mechanism  is  fundamentally  different  in  those  materials  than  in 
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metals.  Comparison  of  the  shock  data  with  more  recent  hydrostatic  data, 
however,  does  indicate  an  appreciable  difference,  so  that  the  elastic- 
plastic  model,  modified  to  include  stress  relaxation,  may  be  at  least  pheno¬ 
menologically  reasonable.  (Fig.  15)  Similar  results  have  also  been  obtained 
for  A1 2O3  in  polycrystall ine  form. 

Whether  these  latter  results  imply  that  brittle  crystals  such  as 
quartz  and  AlgOg  can  deform  extensively  by  dislocation  motion  under  shock  con¬ 
ditions  is  not  clear.  If  so,  the  dislocation  mobilities  must  be  substantially 

If 

I  am  indebted  to  Dr.  E.  B.  Royce  for  pointing  this  out.  R.  A.  Graham  points 
out,  however,  that  the  states  above  the  elastic  limit  for  different  orienta¬ 
tions  of  the  quartz  fall  on  the  same  curve.  Moreover,  there  is  a  large  reduction 
in  electric  polarization  behind  the  second  snock  front.  These  observations  are 
consistent  with  loss  of  rigidity  above  the  elastic  limit. 
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increased  under  shock.  Recovered  samples  from  experiments  of  this  type  are 
usually  in  the  form  of  powder  so  that  brittle  fracture  occurs  somewhere  in 
the  shock  and  rarefaction  process.  Possibly  fracture  is  the  dominant  yield¬ 
ing  mechanism  in  the  compression  part  of  the  wave,  but  intergranular  friction 
maintains  sufficient  shear  stress  to  give  the  stress  differences  observed. 

The  shear  stresses  calculated  for  the  elastic  wave  are  very  much  higher  than 

fil  fi? 

for  metals  and  approach  theoretical  shear  stresses  for  perfect  crystals.  ’ 

It  should  be  noted  that  other  data  still  indicate  a  loss  of 
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rigidity  above  the  yield  point,  notably  in  single  crystal  MgO  and  A^O^. 

Many  materials,  including  metals  and  ionic  and  covalent  crystals, 
exhibit  stress  relaxation.  This  is  most  easily  observed  in  the  decay  of 
the  elastic  precursor  wave  with  distance  of  travel  from  the  impact  surface, 
and  in  the  structure  of  the  wave  between  the  elastic  and  plastic  shock  fronts. 
Stress  relaxation,  or  strain-rate,  effects  are  also  important  in  the  plastic 
shock  front  and,  although  the  structure  of  this  front  can  be  resolved  in  many 
experiments  at  low  stresses,  only  modest  effort  has  thus  far  been  made  to 
correlate  the  shock  structure  with  strain-rate  models. 

From  the  stress  relaxation  relations  (Section  III.D.)  the  plas¬ 
tic  strain-rate  at  the  elastic  front  can  be  determined  from  the  rate  of 
decay  of  the  peak  stress  of  the  elastic  wave.  The  relation  is: 

dP/dt  =  -F/2 

where  F  is  the  plastic  strain-rate  multiplied  by  twice  the  shear  modulus. 

This  can  be  combined  with  the  relation  from  dislocation  theory: 

;  =  (l/2)bNv 

where  y  is  the  plastic  strain  rate,  b  is  the  Burger's  vector,  N  is  the 
mobile  dislocation  density  and  v  the  velocity  of  dislocation  motion. 

Since  the  elastic  front  is  very  steep  it  is  assumed  that  no  time 
is  available  for  multiplication  of  dislocations.  Hence,  N  =  NQ,  the  initial 
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dislocation  density.  If  it  is  further  assumed  that,  following  Gilman,  the 

40 

dislocation  velocities  are  given  by, 

v  =  v°°  exp  -D/t 

where  v”  is  the  limiting  dislocation  velocity,  D  is  a  parameter  called  the 

drag  stress,  and  t  is  the  shear  stress,  the  equations  above  can  be  combined 

and  integrated  to  give  the  peak  elastic  stress  as  a  function  of  distance  of 

fift 

travel  and  impact  stress.  Comparison  of  the  data  then  yields  a  set  of 

compatible  values  for  NQ  and  D  for  a  given  impact  stress. 

The  most  detailed  study  of  elastic  precursor  decay  has  been  per- 
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formed  on  Armco  iron  by  Taylor.  His  results  are  shown  in  Fig.  16.  Reasonable 

values  of  NQ  and  D  do  indeed  give  a  good  fit  to  these  data,  lending  support  to 

the  theoretical  model.  Johnson,  however,  has  recently  pointed  out  that  Taylor 

assumed  that  only  one  slip  system  of  several  possible  systems  was  active  in 

the  (polycrystalline)  iron.^  If  all  systems  are  taken  into  account  through 

an  averaging  process  the  necessary  dislocation  density  is  increased  by  a 

factor  of  about  five.  This  density  seems  somewhat  high  compared  to  that 

obtained  from  independent  measurements;  consequently,  the  validity  of  the 

model  is  still  somewhat  tenuous.  Kelly  and  Gill  is  point  out  that  under  the 

right  conditions  one  might  be  able  to  discriminate  between  various  disloca- 

fift 

tion  models  by  experiments  of  this  type. 

69 

Precursor  decay  has  also  been  studied  in  iron  by  Ivanov,  et  al , 
in  quartzite  by  Johnson, ^  and  in  aluminum  by  Barker,  et  al.^ 

Jones  and  Holland  have  studied  the  effect  of  grain  size  on  the 
Hugoniot  elastic  limit  in  mild  steel. ^  Although  static  tensile  tests  showed 
marked  differences  in  the  upper  and  lower  yield  point  with  varying  grain  size, 
no  effect  on  the  Hugoniot  elastic  limit  was  observed.  They  conclude  that 
under  the  impact-loading  conditions  employed  dislocations  do  not  move  far 
enough  to  encounter  grain  boundaries  --  in  contrast  to  the  case  of  static 
yielding.  The  dynamic  yield  points,  moreover,  were  two  to  three  times  those 
of  the  static  experiments.  They  have  also  observed  Bauschinger  effects  in 
pre-strained  specimens.^ 
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C.  Porous  Solids 


There  is  considerable  interest  in  shock  propagation  in  porous 
solids,  not  only  because  equation  of  state  data  can  be  obtained  over  a  wide 
range  of  densities  and  internal  energies,  as  mentioned  above,  but  also 
because  porous  solids  possess  excellent  shock  buffering  characteristics. 
Hence,  they  can  be  used  for  the  protection  of  structures  from  shock  damage. 

The  collapse  of  pore  space  leads  to  large  losses  of  internal 
energy  as  mechanical  energy.  In  a  steady  state  shock  the  internal  energy 
is  given  by  Eq.  (3): 


E  -  E0  =  (P+P0)(Vq-V)/2  (3) 

The  compressed  specific  volume  is  not  highly  sensitive  to  the  energy;  conse¬ 
quently,  to  a  first  approximation  we  can  neglect  the  energy  dependence  of 
the  P-V  curve  and  visualize  the  energy  loss  as  indicated  in  Fig.  17. 

In  the  solid  material,  with  initial  volume  V$,  a  steady  shock 
to  pressure  P^  carries  the  material  along  the  Rayleigh  line  joining  P^  and 
V$.  The  triangular  area  under  the  Rayleigh  line  represents  the  internal 
energy  of  the  shocked  state.  The  portion  of  this  internal  energy  recoverable 
as  mechanical  energy  is  approximately  the  area  under  the  R-H  curve.  Hence 
the  mechanical  energy  loss  is  the  sliver-shaped  area  between  the  Rayleigh 
line  and  the  associated  R-H  curve.  Clearly,  this  area  increases  substantially 
with  VQ  as  the  porosity  is  increased. 

The  mechanisms  for  energy  loss  cannot  be  precisely  stated  because 
the  solid  is  three-dimensional  on  the  scale  of  the  pore  size.  However,  the 
principal  mechanism  is  probably  initially  the  conversion  of  directed  kinetic 
energy  in  the  propagation  direction  to  acoustic  energy  propagating  in  random 
directions;  various  dissipative  mechanisms  then  convert  this  energy  to  heat. 

Thouvenin  has  proposed  a  one-dimensional  model  (plate-gap  model) 
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for  a  porous  solid  which  allows  no  mechanism  for  energy  losses.  Consequently, 
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he  concludes  that  steady  state  shocks  cannot  exist,  and  predicts  that  a  single 

R-H  curve  will  be  obtained  for  all  initial  densities.  Hofmann,  et  al ,  have 

shown,  however,  that  with  an  energy-dependent  equation  of  state  the  plate-gap 

model  approaches  the  continuum  :teady-state  model  with  increasing  time  and 
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propagation  distance.  The  question  is  therefore  not  which  model  is  correct 
but  under  what  conditions  is  each  model  correct.  Little  experimental  work 
has  yet  been  done  to  establish  the  conditions  under  which  each  model  is  valid. 

The  picture  of  the  process  shown  in  Fig.  17,  although  qualita¬ 
tively  correct,  is  very  much  simplified.  The  R-H  curves  are,  of  course, 
energy  dependent  so  that  each  value  of  VQ  has  its  own  associated  P-V  curve. 

At  lower  pressures  the  collapse  of  pore  space  may  not  be  complete.  Moreover, 
many  porous  solids  exhibit  a  finite  elastic  yield  stress  so  that  a  pressure 
pulse  propagates  as  two  or  more  wavefronts. 

A  more  accurate,  though  still  qualitative,  P-V  diagram  is  shown 
as  Fig.  18.  The  elastic  limit  is  indicated  by  Pg  and  the  pressure  at  which 
the  pore  space  has  completely  closed  is  indicated  by  P  .  The  two  R-H  curves 
above  Pc  are  separated  because  of  greater  shock  heating  of  the  initially 
porous  material.  Partial  compaction  occurs  in  the  region  between  Pg  and  P  , 
and  following  compaction  the  relief  curve  is  substantially  steeper  than  the 
compaction  curve. 
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Linde  and  Schmidt  and  Herrmann  have  proposed  phenomenological 
models  for  the  P-V  relation  of  porous  solids,  including  the  elastic  yield 
point  and  the  compaction  region.  Although  they  require  experimentally  deter¬ 
mined  parameters,  these  models  appear  to  fit  existing  data  reasonably  well. 

The  region  of  partial  compaction  is  least  well  understood,  especially  the 
release  curves  from  a  compressed  state. 

The  pressures  at  which  pore  collapse  is  complete  depend  on  the 

material.  Boade  has  reported  that  compaction  is  complete  in  porous  copper 
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at  pressures  of  about  21  kbar.  In  tungsten  the  compaction  pressure  is 

about  100  kbar  and  in  iron  it  is  26  kbar.  Aluminum  and  graphite  evidently 
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compact  at  lower  pressures,  of  the  order  of  a  few  kbar. 
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Other  models  intermediate  between  the  two  mentioned  have  also 

been  employed.  The  simple  snow-plow  model  assumes  zero  elastic  strength 

(P  )  and  compaction  pressure  (P  ),  and  assumes  further  that  the  compacted 
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material  is  perfectly  rigid.  This  model  is  very  amenable  to  analysis  and 
gives  surprisingly  accurate  results  for  shock  attenuation  in  highly  distended 
materials.  A  modification  which  allows  finite  elastic  strength  but  maintains 
the  assumption  of  perfect  rigidity  of  the  compacted  material  has  also  been 
used  for  approximate  analysis. 

The  elastic  strength  and  the  elastic  wave  speed  are  both  sub¬ 
stantially  reduced  by  porosity.  Butcher  and  Karnes  have  reported  that  the 
elastic  precursor  amplitude  varies  linearly  with  distention  ratio  from  about 
9.5  kbar  in  iron  of  zero  porosity  to  1.5  kbar  at  a  distention  ratio  (i.e. 

ratio  of  initial  density  of  solid  to  initial  density  of  porous  solid)  of 
78 

1.63.  The  elastic  wave  velocity  also  varies  linearly  from  6.0  km/sec  at 
zero  porosity  to  1.15  km/sec  at  a  distention  ratio  of  2.37.  These  results 
seem  to  be  typical  for  most  porous  solids  investigated  to  date. 
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A  third  wave  has  been  observed  in  copper  by  Boade.  The  first 
travels  with  acoustic  velocity  with  an  amplitude  of  a  few  tenths  of  a  kbar; 
the  second  travels  with  a  velocity  of  about  half  the  acoustic  velocity  and 
an  amplitude  of  1.35  kbar.  These  waves  are  followed  by  a  third  that  carries 
the  material  to  the  peak  pressure.  The  origin  of  the  second  wave  is  unclear; 
it  may  be  related  to  a  distinct  stage  of  the  compaction  process. 

The  shock  profiles  observed  in  porous  iron,  graphite,  and  alu¬ 
minum  are  not  steady  but  continue  to  spread  in  the  specimen  thicknesses 

19  78 

employed  in  the  experiments.  ’  The  physical  mechanisms  for  stress  relaxa¬ 
tion  are  not  established;  however,  it  is  possible  that  the  experiments  are 
observing  the  transition  from  the  plate-gap  model  to  the  continuum  steady 
state  model . 


D. 


Plastics 
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Only  a  limited  amount  of  information  is  available  for  plastics 
in  the  stress  range  at  which  one  expects  nonhydrostatic  effects  to  be  appre¬ 
ciable.  Poly  methyl  methacrylate  (Plexiglas),  however,  has  been  studied  by 

several  investigators  and  it  is  reasonable  to  suppose  that  its  behavior  is 
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typical  of  many  polymers.  An  excellent  discussion  of  its  high  pressure 
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behavior  has  been  given  by  Deal. 

There  are  two  distinctive  features  of  the  shock  response  of 

plexiglas  that  are  in  marked  contrast  to  the  behavior  of  metals  or  of  brittle 

solids.  An  elastic  precursor  wave  has  never  been  observed  even  at  low  im- 
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pact  stresses,  and  evidently  one  does  not  exist.  The  observed  shock 

velocities,  moreover,  extrapolate  to  the  longitudinal  sound  speed  at  zero 

stress.  It  is  not  reasonable  to  suppose  that  the  shear  stresses  increase 

indefinitely  with  shock  stress,  although  Ainbinder,  et  al ,  have  observed  an 

increase  by  nearly  a  factor  of  two  in  yield  stress  under  confining  pressures 
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of  about  2  kbar.  It  may  be  that  yielding  occurs  over  a  wide  stress  range 
so  that  a  cusp  in  the  P-V  curve  (Fig.  2)  never  appears.  Measurement  of 
release  curves  or  comparison  of  shock  data  with  hydrostatic  data  would  pro¬ 
vide  valuable  information  about  the  shear  stresses.  Experiments  on  epoxy 
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at  stresses  as  low  as  3.2  kbar  also  show  no  indication  of  an  elastic  wave. 
That  epoxy  cannot  be  treated  as  a  fluid  is  indicated  by  Erkman's  data  on 
shock  decay. (Fig.  19) 

Plexiglas  also  exhibits  pronounced  relaxation  of  the  shocked 

state  with  propagation  distance  in  experiments  in  which  rarefaction  waves 
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from  boundaries  cannot  influence  the  shock  front.  Thus,  it  is  a  highly 
rate-sensitive  material. 

Kolsky  has  devised  a  constitutive  relation  for  polymers  that 

successfully  predicts  the  pulse  shapes  in  plexiglas  and  polyethylene  in  the 
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low  stress  range  in  rod  or  spherical  wave  geometry.  It  is  based  on  a 
viscoelastic  model  and  employs  an  empirically  determined  value,  assumed 
independent  of  frequency,  for  the  phase  angle  between  the  stress  and  strain 
amplitudes  in  a  sinusoidal  wave.  It  would  seem  fruitful  to  attempt  to  apply 
his  analysis  to  plane  waves  as  well.,  and  to  try  to  extend  it  to  include  non¬ 
linear  stress-strain  behavior. 
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V.  CONCLUSIONS 

The  subject  of  one-dimensional-strain  constitutive  relations  applicable 
to  shock  propagation,  with  the  possible  exception  of  high  pressure  equation 
of  state  measurements,  is  clearly  in  its  infancy.  As  experimental  techniques 
become  more  refined  the  concepts  and  assumptions  of  gas  dynamics,  such  as 
discontinuous  shock  fronts  and  equilibrium,  isentropic  rarefactions,  which 
formed  the  early  foundation  for  the  subject  appear  to  lose  pertinency.  Cer¬ 
tainly,  the  assumption  of  isotropic  stress  even  at  stresses  very  much  in 
excess  of  the  yield  stress  has  been  found  to  be  questionable.  Stress  relaxa¬ 
tion  effects  also  seem  to  be  very  common. 

Because  of  the  often  severe  distortion  of  the  shock  structure  by  reflec 
tions  at  interfaces,  the  techniques  which  seem  to  hold  the  best  promise  for 
significant  advances  in  understanding  are  those  that  can  be  imbedded  in  the 
material,  such  as  piezoresistive  gauges,  electromagnetic  velocity  gauges, 
and,  for  transparent  substances,  optical  methods.  Quartz  and  sapphire  gauges 
of  course,  should  be  added  to  this  list  where  the  shock  impedance  of  the 
sample  is  closely  matched.  These,  combined  with  a  general  analytical  treat¬ 
ment  such  as  that  outlined  in  Section  III,  should  provide  a  sounder  base  of 
experimental  information  than  currently  exists. 

The  elastic-plastic  model  has  been  demonstrated  to  give  reasonable  pre¬ 
dictions  for  metals  when  rate  effects  are  not  large  and  at  stresses  only 
moderately  exceeding  the  yield  stress.  Rate  effects  are  the  subject  of  much 
current  research  and  it  may  be  expected  that  progress  will  be  made  in  under¬ 
standing  not  only  such  phenomena  as  the  decay  of  elastic  precursor  waves, 
but,  perhaps  more  importantly,  the  physics  of  the  plastic  shock  transition 
itself.  The  behavior  of  the  yield  stress  at  higher  pressures  is  clearly  not 
well  understood;  whether  or  not  it  increases  substantially  with  pressure  and 
to  what  extent  the  Bauschinger  effect  is  important  are  still  controversial 
topics. 


Brittle  solids  are  even  less  well  understood  and  one  cannot  predict 
with  confidence  whether  for  a  given  material  substantial  shear  stresses  can 
exist  behind  the  "plastic"  shock  front,  nor  what  the  release  paths  from  a 


43 


shocked  state  may  be.  Many  of  these  solids  are  further  complicated  by  phase 
transformations. 

Plastics  appear  to  comprise  another  distinct  class  of  materials  for 
which  a  model  for  the  constitutive  relation  under  shock  needs  formulation. 
Much  more  experimentation  is  required  before  theory  can  fairly  begin  on 
this  problem,  although  the  success  of  Kolsky's  model  at  low  stresses  is 
encouraging. 

Porous  solids  form  still  another  class  in  which  much  work  is  currently 
underway.  The  nature  of  the  compaction  process  is  only  beginning  to  be 
understood  and  is  crucial  to  predictions  of  shock  propagation. 
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TABLE  I 


Stresses  Produced  by  Impact  at  1.5  mm/us.  (kbar) 


Projectile 

Target 

A1 

Fe 

i 

Pb 

W 

Lucite 

A1 

130 

175 

175 

235 

30 

Fe 

175 

280 

275 

390 

70 

Pb 

175 

275 

270 

385 

70 

W 

235 

390 

385 

640 

75 

Lucite 

30 

70 

70 

75 

30 
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Figure  1.  Pressure-Volume  Curves  for  a  Simple  Fluid.  The  Hugoniot  curve 
is  the  locus  of  equilibrium  end  states  attainable  through  a 
series  of  steady  shock  transitions.  The  straight  line  is  the 
Rayleigh  line  for  a  shocked  state,  "C". 
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Figure  2.  Hugoniot  Curve  for  a  Solid  Showing  a  Cusp  at  the  Elastic  Limit. 
A  single  shock  front  is  stable  at  stresses  less  than  "B"  or 
greater  than  "D".  Intermediate  shock  amplitudes  propagate  as 
two  wavefronts. 
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Figure  3.  Stress-Velocity  Curves.  A  shock  in  Material  "A"  with  state 
(P1,u1)  reflects  at  an  interface  with  material  "B"  to  produce 
the  common  state,  (Pg.Ug).  The  curves  A'  and  A  are  the  same 
when  entropy  changes  can  be  neglected. 


-O 

03 

u 

p 

o 

N 

TD 

l/> 

N 

C 

-Q 

3 

fO 

fO 

E 

O 

0) 

O 

•M 

x: 

o 

•M 

C 

>> 

4- 

03 

"ol 

O 

p- 

03 

N 

c 

p 

N 

*4- 

p— 

3 

*0 

E 

l/> 

03 

03 

T3 

03 

x: 

•p 

x: 

4-> 

r  " 

4-> 

l/i 

O 

O 

4-> 

C 

4-> 

=3 

03 

CT> 

>0 

cn 

S 

a; 

TD 

E 

u 

*r 

•r 

C7> 

“D 

•M 

•p 

C 

U 

00 

03 

TD 

C 

03 

03 

03 

x: 

■M 

> 

• 

h- 

C 

03 

03 

3 

U 

00 

Q 

Q. 

•p 

E 

o 

r* 

C 

* 

00 

C 

C 

3 

•p 

2 

00 

03 

O 

4-> 

u 

03 

U 

3 

t/> 

o> 

03 

<o 

U 

x: 

o 

o 

U 

4-> 

4-> 

=) 

>0 

oo 

03 

U 

2 

x: 

p 

h- 

03 

fO 

4- 

JZ 

•p 

o 

h- 

4-> 

• 

c 

x: 

u 

•r* 

, 

2 

4-> 

O 

E 

00 

03 

O 

O 

X3 

03 

O 

3 

00 

U 

C 

00 

2! 

=3 

o> 


i 


j 


CUR  RE 
LEADS 


Figure  5.  Manganln  Pressure  Transducer.-  Four  pi-shaped  manganin  wires  are 
shown  Imbedded  In  epoxy.  The  change  in  voltage  Is  monitored  as 
the  shock  passes  over  the  wires.  Current  is  held  constant. 
(After  Keough^) 


Figure  7.  Schematic  of  Electromagnetic  Velocity  Technique.  Motion  of  foil 
imbedded  in  insulator  is  monitored  by  observing  EMF  developed. 
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Figure  8.  Laser  Interferometer.  Fringes  developed  by  interference  of 
beams  reflected  from  moving  and  stationery  mirrors  record 
displacement  of  target  free-surface.  (After  Barker*^) 


TARGET 


Figure  9.  Velocity  Interferometer.  Fringes  are  developed  by  interference 

of  reflected  beam  with  that  reflected  from  same  surface  earlier 

in  time.  Fringe  counting  rate  is  proportional  to  acceleration 
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of  surface.  (After  Barker  ) 
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Crystal  (001  orientation).  Thickness  2.133  run.  Impact  with  A1 
projectile  produced  step  stress  pulse  of  25  kbar  in  LiF.  (After 


Figure  12.  Stress-Volume  Curves  for  One-Dimensional  Strain  Based  on  Elas¬ 
tic-Plastic  Model.  Relief  of  stress  from  a  (shocked)  state  "C" 
follows  path  CDE.  PQ  Is  Hugonlot  elastic  limit,  Y  Is  yield 
stress. 
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Figure  13.  Comparison  of  Hugoniot  and  Hydrostatic  Compression  States  in 
2024  Aluminum.^ 
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Figure  14.  Yield  Stress  in  2024  Aluminum  Shocked  to  130  kbar.  Values  are 
determined  from  comparison  of  compression  and  relief  states. 
(After  van  Thiel  and  Kusubov^®) 


Figure  15.  Comparison  of  Hydrostatic  and  Shocked  States  in  a-quartz 

61  -65 

(x-cut).  "  Disparities  in  the  hydrostatic  curves  cause 
uncertainties  in  the  shear  stress  of  the  shocked  state  above 
80  kbar. 
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Figure  16.  Decay  of  Precursor  Amplitude  in  Armco  Iron.  Curves  are  fitted 
using  Gilman  dislocation  velocity  model.  (After  Taylor^) 


Figure  17.  Approximate  Hugoniot  and  Rayleigh  Lines  for  Solid  and  Porous 
Material.  Cross-hatched  areas  represent  energy  dissipated  in 
shock  transition. 
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Figure  18.  Comparison  of  R-H  Curve  of  Porous  Material  with  Hydrostat  of 
Solid  Showing  Hugoniot  Elastic  Limit,  Compaction  Region  and 
Stress  Difference  due  to  Shock  Heating. 


Figure  19.  Decay  of  Peak  Shock  Stress  with  Propagation  Distance  in  Epoxy. 

Solid  curve  is  predicted  on  the  basis  of  fluid  model  neglecting 
shear  stresses.  X  is  the  thickness  of  the  flyer  plate  that 

0  cp 

produced  the  shock.  (After  Erkman  ) 
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